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of the reduced density matrix of the two-qubit system after integrating out the elec- 
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lowest energy states, and Class B states which are the linear combinations of the 
symmetric and the antisymmetric Bell states. Our results for an initial Bell state are 
similar to those obtained before for the same model derived under the Born-Markov 
approximation. However, in the Class A states the behavior is qualitatively different: 
under the non-Markovian evolution we do not see sudden death of quantum entan- 
glement and subsequent revivals, except when the qubits are sufficiently far apart. 
We provide explanations for such differences of behavior both between these two 
classes of states and between the predictions from the Markov and non-Markovian 
dynamics. We also study the decoherence of this two-qubit system. 
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I. INTRODUCTION 

Quantum decoherence [l| and entanglement between one two-level atom (2LA) and 
an electromagnetic field (EMF) has been treated by many authors (see, e.g., the review j^) 
including the present ones 0, o, 6]. Zyczkowski et al j3] have studied discrete dynamics of 
entanglement of a bipartite quantum state. For two 2LA-EMF systems, each of the two 
atoms can be assumed to interact only with its own cavity EMF, or with a common EMF 
(of course they can also interact with each other). The case of each atom interacting with 

riQ 

its own field was studied by Yu and Eberly [8|, 19| who reported on the appearance of 'sudden 
death', or finite time disentanglement. The case of two atoms interacting through a common 
field was studied by Ficek and Tanas under the Born-Markov approximation (BMA) |10i]. In 
contrast to the 'sudden death' found in two qubits in disjoint EMFs they show the existence 
of dark periods and revival of quantum entanglement. The present authors have studied 
the case of two 2LA interacting through a common electromagnetic field without invoking 
the BMA HI]. This essay is a synopsis of this unpublished work, focussing on the new 
non-Markovian behavior of entanglement dynamics and decoherence features. 

nn 

In our prior work [4, |5| , we used the influence functional formalism with a Grassmannian 
algebra for the qubits (system) and a coherent state path integral representation for the 
EMF (environment). Here an operator method is used in conjunction with perturbation 
theory. Because the assumption of an initial vacuum state for the EMF allows a full re- 
summation of the perturbative series, a closed expression for the evolution of the reduced 
density matrix of the two qubits can be obtained. This approach incorporates the back- 
action of the environment (field) on the system (two 2LA) self-consistently which (almost 
always) engenders non-Markovian system dynamics the behavior of which cannot be fully 
comprehended if one imposes the restrictions of BMA. 

In this broader context we can understand certain effects like 'sudden death' Sj as conse- 
quences of rather special arrangements: Each atom interacting with its own EMF precludes 
or lessens the chance that the atoms may be entangled through the mediation of the EMF. 
The separate field case studied by Yu and Eberly corresponds to the limit of the common 
field case studied here when the qubits are separated by distances much greater than the 
correlation length characterizing the total system. Our study shows that for a wide range 
of spatial separations within the correlation length, entanglement is robust and there is 
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lOj highlights the 



no sudden death. The difference between our results and that of Ref. 
non-Markovian effects in the evolution of quantum entanglement. 

Under the usual two level, dipole and rotating wave approximations (RWA) but with- 
out invoking the Born-Markov approximation (BMA) we derive a non-Markovian master 
equation which has a more elaborate structure than the usual Lindblad form: it contains 
extra terms that correspond to off-diagonal elements of the density matrix propagator. We 
observe very different behavior in two classes of states, superpositions of highest and lowest 
energy states which we call Class A states, and the usual antisymmetric |— ) and symmetric 
1+) Bell states [12;i], which we call Class B states. We find similar behavior in the Class 
B (Bell) states compared with the findings of Ref lOj but qualitatively different behavior 



in the evolution of Class A states. Ficek and Tannis [10| found that their evolution leads 
generically to sudden death of entanglement and a subsequent revival. In our more complete 
treatment of the atom-field dynamics we see the former effect present for large values of the 
inter-qubit distances. However, sudden death is absent for short distances, while there is no 
regime in which a revival of entanglement can take place. (Both our results and the ones of 
Ref. jlO| rely on the Rotating Wave Approximation, which is not adequate for very large 
interqubit distances - details in sections below.) 

This paper is organized as follows: Section 2 contains the main derivation. We write 
down the Hamiltonian for two 2-level atoms interacting with a common electromagnetic 
field at zero temperature, and we compute the relevant matrix elements for the propagator 
of the total system by resummation of the perturbative series. We then determine the 
evolution of the reduced density matrix of the atoms, which is expressed in terms of seven 
functions of time. We compute these functions using an approximation that amounts to 
keeping the contribution of the lowest loop order for the exchange of photons between 
the qubits. In Section 3 we examine the evolution of the reduced density matrix for two 
classes of initial states. We then describe the time evolution of quantum entanglement in 

n 

these states via the concurrence function [13] with spatial separation dependence plotted 
for some representative cases. We discuss the results on disentanglement, describe the non- 
Markovian dynamics and identify the source of difference from quantum dynamics under the 
Born-Markov approximation. In Section 4, we study the decoherence of this system when 
the two qubits are initially disentangled. 

In the present synopsis we shall describe our approach, the main steps in the deriva- 
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tions and the physical significance of our findings but leave many details out as they can be 
found in [ll]. As a literature update (from the time ofJll|) on nonMarkovain entanglement 
dynamics we mention the work of the Turku group iJ] on atom-field interaction with ex- 
perimental verification possibilities and that of Paz and Roncaglia 15|, Lin and Hu [l^ on 
two harmonic oscillators (detectors) interacting with a common EMF. The latter paper ad- 
dresses a wider scope of issues including outside-light cone entanglement, and entanglement 
generation from initial separable states with separation dependence. Lin and Hu find that 
when two initially entangled detectors are still outside each other's lightcone, the propor- 
tionality of the degree of entanglement to the spatial separation oscillates in time. When the 
two detectors begin to have causal contact, an interference pattern of the relative degree of 
entanglement (compared to those at spatial infinity) develops a parametric dependence on 
d. The detectors separated at those d with stronger relative degree of entanglement enjoy 
longer disentanglement times. In the cases with weak coupling and large separation, the 
detectors always disentangle at late times. For sufficiently small d, the two detectors can 
have residual entanglement even if they initially were in a separable state, while for d a little 
larger, there could be transient entanglement created by mutual influences. 



II. TWO 2-LEVEL ATOMS INTERACTING VIA A COMMON 
ELECTROMAGNETIC FIELD 

A. The Hamiltonian 

We consider two 2- level atoms (2LA), acting as two qubits, labeled 1 and 2, and an 
electromagnetic field (EMF) described by the free Hamiltonian Hq 

Ho = hJ2uj,blK + huoS^^^s^^^^ + hooj^^^s^y (1) 

a 

where a = (k, a) is a collective index running over momenta k and polarizations a = 1,2 
of the photon bath, oOa = t^k is the frequency of the k electromagnetic field mode and 
uJo the atomic frequency between the two levels of the atom, assumed to be the same for 
the two atoms. The electromagnetic field creation (annihilation) operator is 6^ (ba), while 
(S'i"'') are the spin raising (lowering) operators for the atom. We will define the 
pointing vector from 1 to 2 as r = r2 — ri and we will assume without loss of generality that 
ri + ra = 0. 
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The two 2LAs do not interact with each other directly but only through the common 
electromagnetic field via the interaction Hamiltonian 

a 

While the coupling functions Qa = g-ka can be completely general, its relevant form for 
the dipole coupling of the electromagnetic field to the dipole of the atom d is 

ga ■= g^a = ^==d ■ ek^, (3) 

where ekcr are the polarization vectors, Ad is the transition matrix element of the momentum 
p of the atomic electron (assumed reals) , d being the corresponding unit vector; V is the 
volume of space. We have assumed that Ad is the same in both atoms. The total Hamiltonian 
of the atom-field system is 

H = Ho + Hi. (4) 

The atom-field interaction Hamiltonian ([2]) is derived under the dipole and the rotating 
wave approximations (RWA) (see Appendix A in Ref. ^]). The RWA keeps only the terms 
in the interaction-picture Hamiltonian that correspond to resonant coupling but ignores all 
the rapidly oscillating terms. For a single qubit system the RWA is self-consistent. However, 
in the two-qubit system since we keep terms in the Hamiltonian that vary in space as e^^'^ 
the RWA is consistent only if r << t. This condition is satisfied in most realistic situations, 
but it fails in the formal limit r —>■ 00. 

The two level approximation holds only for photon frequencies that are not much larger 
than Uo. For this we introduce a high-frequency cut-off such that cUoC << 1 [4]- In 
what follows, the cut-off is implemented by inserting a factor e"'^'^'' in all integrations over 
momenta k that appear in the definition of the evolution functions. For t,r » e^^, the 
cut-off dependence can be absorbed in a renormalization of the basic parameters. 

As the results in Section II B to II D do not depend on the explicit form of the coupling 
constant Qa (and they can also be applied to non-optical systems with similar coupling) we 
leave the index a completely general and only use the specific form ([3]) in Sec. II E. 



B. Perturbative expansion and resummation 



We assume an initial factorized state at t = for the combined system of atoms+field 
of the form \0) (S> {ip), where \0) is the vacuum state of the EMF and lip) is a vector on 
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the Hilbert space of the two 2LA's. To study the dynamics of this two 2LA-EMF system 
we seek the action of the evolution operator e~*^* on such vectors. The calculations are 
simplified by using the resolvent expansion of the Hamiltonian 

27r J E-H + if) ^ ^ 

which allows for a perturbation expansion 

{E - H)-' = {E- Ho)-' + {E- Ho)-'Hi{E - Ho)-' 

+ {E - Ho)-'Hi{E - Ho)-'Hi{E - Ho)'' + ... (6) 

Of relevance for the computation of the reduced density matrix of the two qubits are 
matrix elements of the form {z; i',j'\{E — H)^^\0; i,j), where z refers to a coherent state of 
the EM field and i,j take on the values of 0, 1 corresponding to the ground state of a single 
qubit and to the excited state respectively. We compute the matrix elements above through 
the perturbation expansion ([6]). A resummation of the perturbative series is possible, and 
it leads to an exact expression for the matrix elements of the resolvent-see Appendix A 
of Ref. jll| for details. An inverse Fourier transform yields the following values for the 
non- vanishing matrix elements of the propagator. 



a 

dEe-'^^ 



(2;0, l|e-^-^*|O;0, 1) 



1 



E-ujo-a{E) - (3{E,r) 



+ 



1 



(;2;l,l|e-^^*|O;0,l) 

(2;0,0|e-^-^*|O;0,0) 
(^;l,l|e-^^*|0;l,l) 

(;z;0,0|e-^-f^*|O;l,l) 

(2;0,l|e-^^*|O;l,l) 
(2;l,0|e-^^*|O;l,l) 



E~uJo-a{E)+ f3{E,r) 



-■■ v+{t) 



dEe 



-iEt 



E-iUo-a{E) - l3iE,r) 



E-uJo-a{E)+(3{E,r) 



-■■ v-{t) 



dEe 



iEt 



E-2uJo-2a{E-uJo)- f{E,r) 



--: u{t) 



dEe~^^* 

^E-2uJo~ 2a{E - Uo) - f{E, r) 



e 2 



ya{t) 



\ e'-Kit) 



(7) 



(8) 



(9) 
(10) 
(11) 

(12) 
(13) 
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where, 



Sa{t) 


-J- 






[ Kit) ) 




a{E) : 


a 


PiE,r): 


= E 

a 



The matrix L is 



where 



1 



dEe 



-iEt 



(E-LUo- aiE) - PiE, r)e^-){E - uj^) 



dEe 



iEt 



I 



E - 2uJo 



E(i-^) 



ab 



at 



\ E—u)o—uib 



9l 



9a g*''''" 



E -Ur, 



L :-- 



E e 
e s 



e 



E — Uo — uj, 
1 



ab 



a{E - UJa)5ab + gagb' 
{(3{E - UJa,r)6ab + QaQbi 



+ 



gj(k-k')r 



E — 2uJo E — Ua — UJb 

1 1 

+ 



E — loJo E — Ua — UJh' 



(14) 

(15) 

(16) 
(17) 



(18) 



(19) 
(20) 



E -Uo-UJa V 

and the overbar denotes complex conjugation. 

The exphcit definitions of the kernel Hab and of the function / can be found in [ll|] but 
will not be needed for the issues we are interested in here. 



C. The reduced density matrix 

We next compute the elements of the reduced density matrix for the qubit system by 
integrating out the EM field degrees of freedom 

P^f(t)= E f^Sa^) /HM^*](0;aie^'^*k;^',/)(^,^,j|e-^'^*|0,^o,Jo), (21) 

io,jo,i'od'o 

where [dz] is the standard Gaussian integration measure for the coherent states of the EM 
field. 

Substituting Eqs. flTl fT^ into fl^T]) we obtain through a tedious but straightforward cal- 
culation the elements of the reduced density matrix 
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where 



p\{t) = p^m^m (22) 

Pllit) = plliO)uit)vlit) + pll{0)u{t)vlit) (23) 

pIoW = plliOMt)vlit) + plliO)uit)v*_it) (24) 

Poo(^) = plMt) (25) 

PooW = Poo(0)^+(t) + Poo(0)t;-(i) + Poi(0)pi(t) + Plli0)p2it) (26) 

pM = Pl»Mt) + Pl»v-it) + Pll{0)pm + p}j(0)/il(t) (27) 

P'olit) = pS}(0)|t;+r(t) + P?J(0)t;4tKW + p}^(0)|t;-P(t) 

+ pj?(0)t;_(t)<(t) + p}l(0)/ti(t) (28) 

P?oW = P?o(0)|^^+r(t)+P^?(0)|t;_r(t)+p°}(0)t;4t)t;*(t) 

+ pl',iO)v4t)vlit) + plliO)K,it) (29) 

pSSW = l-PnW-P^iW-pl^W (30) 



p,{t) = Y.9aVa{t)sl{t) (31) 

a 

P2{t) = Y.9aVa{t)sl{t)e-''^- (32) 

a 
a 

«:2(t) = E^a(t)i^;*(t)e-^'^-, (34) 

a 

and the functions u{t),v±{t) were defined in Eqs. (fTT]) . (jS]) and ([21) • 



D. Explicit forms for the evolution functions 

Eqs. fl22ti30p provide an exact expression for the evolution of the reduced density matrix 
for the system of two qubits interacting with an EM field in the vacuum state. The evolution 
is determined by seven functions of time m, v±^ ki,2, Pi,2 defined above. To obtain the explicit 
forms for these functions involves summing over all modes labeled by a, namely, performing 
an integration over all momenta k and summing over all polarizations. 

We employ the following approximations: 
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1. We assume weak coupling (A^ << 1) and ignore the contribution of all processes that 
involve the exchange of two or more photons between the two qubits. 

2. We ignore branch cut terms that appear from the inverse Fourier transform of the 
coefficients and keep only the contribution of the dominant poles. One may prove that 
the branch-cut term only becomes significant at the very- long-time limit, well beyond 
the time scales associated to decay. 



1. The functions u, v± 

With the approximation above, the contribution of the function / drops out from the 
definition of u. Thus we obtain 



u 



E-2iUo- 2a{E - Uo) 



v± 



dEe 



i.Et 



1 



± 



1 



E-iUo- a{E) - (3{E, r) E - lUo - a{E) + /3{E, r) 



(35) 
(36) 



For small A there is only one single pole and thus they become 



p-iujot-Vat 

^ f-iat-Vrt _j_ ^iat+Trt 



(37) 
(38) 



In the equations above, we renormalized the frequency uJq by a constant divergent term 
that arises from a{E) (see j4|). The parameters ro,rr and a{r) are defined as 



Fo := —Ima{ijJo) 
-a{r) + iVr := P{u}o,r). 



(39) 
(40) 



For the coupling 



and Eqs. (fT6] - [T7|) we find 



To 
F. 



2nr 



(sin LU^r + 



cos UJoT Sm UoT , 



UJoT 



2^2 



UJqV 



Id ■ r) [smuJor + 6 

uJnr 



, sm UpT , 



(41) 



(42) 



where f is the unit vector in the direction of r. 
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The term a{r) is a frequency shift caused by the vacuum fluctuations. It breaks the 
degeneracy of the two-qubit system and generates an effective dipole couphng between the 
qubits. The constant Tq corresponds to the rate of emission from individual qubits; it is the 
same as that from a single qubit interacting with the electromagnetic field. The function 
Tr is specific to the two-qubit system. It arises from Feynman diagrams that involve an 
exchange of photons between the qubits. As r — > 0, — Fq and as r — > oo, F^ 0. Note 
that the ratio F^/Fq, while smaller than unity, is of the order of unity as long as r is not 
much larger than uj~^. 



2. The functions Ki^2{'t) 



We first compute the functions z/^, v'^ of Eq. f|T5|) keeping terms up to second loop order 



Kit) 



dEe 



-lEt 



E 



9b 



6. 



ab -|- ^ab + ©afc 



where S and are given by Eqs. f|T9l) and fl20|) . To leading order in we obtain 

dEe-''^' 

Mt) = I'aW = 9a 



(43) 



(•44) 



[E - 2uJo - 2a{E - cu„)] [E - - uj^ - a{E - Ua) - I3{E - r))| 
Keeping only the pole terms from the integral in (jUj) we obtain 



Va{t) = 9a 



(45) 



— cUfl — cr — iTq + iV 

We then substitute the expression above for z/^ into Eqs. ( |33l) and ( IMll . Using the values of 
Eq. ([3]) for Qa-, and summing over the polarizations and angular variables, we obtain 

A2 



«i(t) 

«2(t) 



-2Tat 



37r2 
A^ 



g"""^"' + g"'^''* - 2e-(^"+^'-)* cos[(cuo - fc - (y)t\ 
(/t - a;. + a)2 + (Fo - F.)2 ^^^^ 



-2rot l"^ + e"^^'-* - 2e"(^o+^'-)* cos[(u;o - A; - a)t\ 







(k-oOo + ay + iTo-Tr)' 



(sin kr + 



cos kr sin A;r , 



(d ■ f)^(sin /cr + 3 



cos kr sin /cr , 



kr 



(47) 



To compute these functions we use an approximation scheme of replacing the Lorentzian 
in fl46tl47p by a delta function, which is valid for time-scales t » u;~^, i.e.. 



TT 



{k-uJo + af + (Fo - Trf Fo - F 



-5{k - uJo + a) + O 



Fn - F. 



0J„ 



21 



(48) 
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To leading order in we obtain the simple expressions 



K,{t) = ToK{t) (49) 

K2{t) = TrK{t) (50) 



where 



^(t) ^ ,^_^e-2ro*(e-rot _ ^-r^ty^ (51) 

III. DISENTANGLEMENT OF TWO QUBITS 

In this section, we employ the results obtained above to study the evolution of the two 
qubits initially in an entangled state. We shall focus on the process of disentanglement 
induced by their interaction with the field. 

A. Class A states: Initial superposition of |00) and |11) 
We first examine the class of initial states we call Class A of the following type 



\^o) = ^Jl-pm + y^\ll), (52) 

where < p < 1. Recall our definition |/) = |11) and \0) = |00). From Eqs. ([22]- [30]) we 
obtain 



-pMt) - K,it)]\-){-\ + K,{t)]\+){+\ + [1 -pV^r»* - 2p^,mO){0\, (53) 



where the functions Ki(t) and K2{t) are given by Eqs. (1191450 



These results are quite different from those reported in Ref. [10|, which were obtained 
under the Born-Markov approximation. While the |/)(/| and |/)(0| terms are essentially 
the same, the |— )( — | and |+)(+| ones are not, as they involve non-diagonal elements of the 
density matrix propagator. For comparison, we reproduce here the explicit form of these 
matrix elements in our calculation 

pX = Pf^f""''(^-"'' - ^-"''f (54) 
p-_ = pe-2ro*(e-ro*-e-r'-*)2, (55) 
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and in that of Ref. lOj (translated into our notation): 



plit) = p^^e-^^«*(e"^^^* - e-^^«*) (56) 

i — i r 

p--it) = p^^e-'^^\e'^^' - e-2^0*). (57) 

J- + i r 

For large values of r, << Fq the expressions above coincide. However, for smaller values 
of r their difference is significant. We note that even though r > as the regime ujor << 1 
is physically meaningful, as long as UodB << 1. In this regime the difference between ( !55ll 
and ( 1571) is substantial: Fq — F^ and at times Fot ~ 1 we obtain (Fq — Tr)t « 1. According 
to the Markovian results of in this regime the |+)(+| term is of order 0(A°) and hence 
comparable in size to the other terms appearing in the evolution of the density matrix. 
However, according to our results, which are based on the full non-Markovian dynamics, the 
|+)(+| term is of order ^-^^ and hence much smaller. 

In general, for ujQr << 1, we find that the |— )(— | and |+)(+| terms contribute little to 
the evolution of the reduced density matrix and they can be ignored. Since these terms are 



1Q|, we 



responsible for the sudden death and subsequent revival of entanglement studied in 
conclude that these effects are absent in this regime. 

B. Bell states 

For the case of an initial Bell state [12] |±) for the two qubits we find 



We see that the symmetric |+) state decays super-radiantly with decay rate Fq + F^ and 
the anti-symmetric state |— ) decays sub-radiantly with the rate Fq — F^. The results agree 



qualitatively with those obtained in Ref. lOj under the Born-Markov approximation. 



Non-Markovian Features 



Let us try to understand these results with some discussions. 
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1. Differences from Born-Markov 

For an initial Class A state (152|) . the |/) component decays to the vacuum, but it also 
evolves into a linear combination of antisymmetric |— ) and symmetric |+) Bell states. How- 
ever, if the qubits are close together the evolution to Bell states is suppressed. This be- 
havior is qualitatively different from that of Ref. lOj, which was obtained through the 
Born-Markov approximation. The corresponding terms differ substantially at the physically 
relevant time-scales. As a consequence, we find that there is neither sudden death nor revival 
of entanglement in this regime. 

In order to explain this difference we note that the Born-Markov method involves two 
approximations: i) that the back-action of the field on the atoms is negligible and ii) that all 
memory effects in the system are insignificant. When the qubits are found within a distance 
much smaller than their characteristic wavelengths, it is not possible to ignore back-action. 
The virtual photons exchanged by the qubits (at a rate given by r^.) substantially alter the 
state of the electromagnetic field. 

On the other hand, the effect of virtual photons exchange between qubits drops off quickly 
at large separations r - the two qubits decay almost independently one of the other. Hence, 
the Born-Markov approximation - reliable for the case of two separate qubits each interacting 
with an individual field - also gives reasonable results for the two qubits interacting with 
a common field. In this regime sudden death is possible, but not revival of entanglement. 

w 

In this sense, our results effectively reduces to those of Ref. [.8!]: when the distance between 
the qubits is much larger than any characteristic correlation length scale of the system it 
behaves as if the two qubits are found in different reservoirs. 



2. The origin of the non-Markovian behavior 

From Eqs. (!22] - [30l) and ( |37l 1381) we note that the diagonal terms of the reduced density 
matrix propagator all decay exponentially, which is a characteristic sign of Markovian behav- 
ior. Hence, as far as this part of the evolution is concerned, our results are fully compatible 
with the Markovian predictions. 

However, the behavior of the non-diagonal terms in the reduced density matrix propagator 
is different. Eqs. fl22ti30p show that the only non-zero such terms are ones that describe 
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the effect of successive decays, for example the |11) state first decaying into |— ) and then 
|— ) decaying into the ground state |00). Hence, the pl(t) term consists of one component 
that contains the remaining of the | — )( — | part of the initial state and another component 
that incorporates the decay of the part of the initial state towards the state |— ). In 

our calculation, the latter term is encoded into the functions Ki^2{t), which are obtained by 
squaring the amplitudes z/k(t) as in Eqs. (I3314MI) . The amplitudes z/k(t) are obtained from 
the summation of two Feynman diagrams -see Eq. (jSj). The structure of the poles in Eq. 
(144|) reveals that the first Feynman diagram describes the decays of the |11) state, while the 
second one corresponds to processes involving the |01) and |10) states. 

The key point is that the evolution functions ki,2(^) ai'e squares of the corresponding 
amplitudes, and for this reason they contain interference terms between the two types of 
processes that contribute to z^k(^)- Major difference of our non-Markovian results stems from 
the constituents of the off-diagonal terms. In the Markov approximation these term involve 
summation (subtraction) of probabilities rather than of amplitudes; hence, it ignores the 
interference between these two processes. 

IV. DECOHERENCE OF TWO QUBITS 

We now turn our attention to the decoherence issue. We want to compare the rate of 
decoherence of the two qubit system by the EMF with that of one qubit alone in order to 
understand how the decoherence of one qubit is affected by the presence of another. We 
assume that the initial state of the two 2LA-EMF system is of the form 

(vp|1) + /l^lO)) ® |0) = v^llO) + /r^lOO). (59) 

where the first qubit is prepared initially in a superposition of the |0) and |1) states, and the 
second qubit lies on the ground state. From Eqs. ( [22] - [30l) . we obtain the density matrix of 
the combined qubit system 

p(t) =p(|t;+niO)(10| + |t;_n01)(01| +t;*t;+|01)(10| +t;_t;;|10)(01|) 
+ ^Jp{l~p) (<|10)(00| +t;l|01)(00| +t;+|00)(10| +t;_|00)(01|) 

+ (l-p(|t;+|2 + |t;_n)|00)(00|, (60) 

where the functions v±{t) are given by Eq. (!38|) . 
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The two qubits become entangled through their interaction via the EM field. To study 
the decoherence in the first qubit, we trace out the degrees of freedom of the second one, 
thus constructing the reduced density matrix 

Ut) = + ^jpii-p{v+\o){i\ + vi\m) + (i -p\v+\') io)(o|. (61) 

At large interqubit separations = = (T(r), whence f+ ~ ^-'^'^ot-rot ^ ^^le off-diagonal 
elements decay within a characteristic time-scale of order Fq"^. These results coincide with 
those for the single qubit case-see Refs. . However, for smaller values of r, the results are 
substantially different. The entanglement with the second qubit leads to a departure from 
pure exponential decay. In particular, for uJor << 1, F.^. ~ Fq. This implies for times longer 
than Fq ^ a substantial fraction of the off-diagonal elements persists. This decays eventually 
to zero within a time-scale of order [Fq — F,.]^^ >> Fg"^. Hence, the qubit preserves its 
coherence longer when there is another quibit close by. 

The reduced density matrix of the second qubit is 

kit) = p\v.\'\l){l\ + ^p{l-p) {v-\0){l\ + v*_\l)0\) + (l -p|i;-r) |0)(0|. (62) 

Note that at small inter-qubit separations the asymptotic behavior ( for Fot >> 1) of Pi{t) 
is identical to that of P2{t)- The second qubit (even though initially in its ground state) 
develops a persistent quantum coherence as a result of the interaction with the first one. 

It is also of interest to consider thermal field environment. Shresta et al J5| have 
studied the decoherence of one qubit in a finite temperature EMF. Yu and Eberly [l7| have 
considered the separability of the joint state of a collection of two-level systems at finite 
temperature. They conclude that since only separable states are found in the neighborhood 
of their thermal equilibrium state unimpeded thermal decoherence will destroy any initially 
arranged entanglement in a finite time. 
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